An A r threesome: 
Matrix models, 2d CFTs and 4d M = 2 gauge theories 

Ricardo Schiappa" and Niclas Wyllard 6 

a CAMGSD, Departamento de Matematica 
Instituto Superior Tecnico 
Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal 

b Department of Fundamental Physics 
Chalmers University of Technology 
S-412 96 Goteborg, Sweden 

(N ; schiappa@math.ist.utl.pt, n.wyllard@gmail.com 

Abstract 

We explore the connections between three classes of theories: A r quiver matrix 
models, d = 2 conformal A r Toda field theories and d — 4 N — 2 supersymmetric 
conformal A r quiver gauge theories. In particular, we analyse the quiver matrix 
models recently introduced by Dijkgraaf and Vafa [T\ and make detailed comparisons 



> 



with the corresponding quantities in the Toda field theories and the M = 2 quiver 
gauge theories. We also make a speculative proposal for how the matrix models 
should be modified in order for them to reproduce the instanton partition functions 
in quiver gauge theories in five dimensions. 



1 Introduction 



The AGT relation [2J, which is a relation between Nekrasov partition functions 
[3] in (conformal) d = 4 J\[ = 2 quiver gauge theories and correlation functions 
in conformal field theories in two dimensions, has been studied in several papers 
over the past couple of months. The original conjecture [2] involves a relation 
between Nekrasov partition functions in d = 4 SU(2) (or Ai) quiver gauge theories 
[1] and correlation functions in the Liouville d = 2 conformal field theory. It was 
subsequently extended [5j to a relation between the A r quiver gauge theories [1] and 
the d = 2 conformal A r Toda field theories. The proposals in [21 [5] have passed many 
non-trivial checks, see e.g. [SI El IE]- In [9] the A\ AGT relation was extended by the 
inclusion of surface, Wilson and 't Hooft operators in the A\ quiver gauge theories 
and proposals were made for the corresponding quantities in the Liouville theory. 
Another line of investigation concerns the extension to non-conformal SU(2) gauge 
theories [10]. A suggestion for how this approach should be modified to capture 
the instanton partition function for pure SU(2) gauge theory in five dimensions 
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was presented in [TT]. In [12], another relation between four-dimensional and two- 
dimensional theories was uncovered which is similar in spirit to the AGT relation. 

Recently Dijkgraaf and Vafa presented an argument explaining the AGT relations 
PQ0. The argument involves relating the relevant quantities in the two theories via 
an intermediate matrix model. The first step is to realize the gauge theories in 
string theory using geometric engineering [15] and then use the relation to matrix 
models via a large N duality, together with the relation between matrix models 
and conformal field theories [16j E] to recover the AGT relation. This chain of 
arguments removes some of the mystery of the AGT relation. More importantly, it 
also implies that there are now three different ways to compute the same quantities: 
using the 4d quiver gauge theories, using the 2d Toda theories, or using the Qd 
quiver matrix models. In all three cases a Riemann surface plays a crucial role: in 
the gauge theory the Riemann surface is related to the Seiberg-Witten curve, in the 
Toda theory the Riemnan surface is the manifold on which the theory is defined, 
and in the matrix model the Riemann surface is the spectral curve arising from the 
loop equations in the large N limit. 

The goal of this paper is to develop and exemplify how calculations are performed 
in the matrix model framework. We will rederive several known results in quiver 
gauge theories and Toda field theories from the matrix model integrals. We also make 
a speculative proposal for how the matrix models should be modified to reproduce 
the Nekrasov partition function for quiver gauge theories in five dimensions. 

In the next section we review the AGT relation for the case of the A r theories and 
in section [3] we describe the A r quiver matrix models introduced in [TJ . In the two 
subsequent sections we then perform several matrix model calculations and compare 
the results with the Toda theories and the quiver gauge theories: In section H] we 
treat the Ax model and in section we discuss the A r models for general r. Finally, 
in section Owe describe our proposal for how the matrix models should be modified 
in order to describe the Nekrasov partition function for quiver gauge theories in five 
dimensions. In the appendix some technical details are collected. 

Note added: After this paper was finished [HI [19] appeared which have some 
overlap with some parts of this paper. 

2 The A r AGT relation 

In this section we review the AGT proposal for the class of theories based on the 
A r Lie algebras. We start with a brief recap of the A r Toda field theories, followed 
by a summary of the A r quiver gauge theories and then describe the AGT relation 
connecting the two classes of theories. 

1 Another argument, using M-theory, for the validity of the AGT relation was presented in |13j : 
see also the follow-up paper [14] . 
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2.1 The A r Toda field theories 

The A r Toda field theories are defined by the action 

JV-l 



S= [ dVVs ^9 ad (da^ 9,0) + p. £ e b <^> + 



(2-1) 



where g a d (a,d = 1,2) is the metric on the two-dimensional worldsheet, and R is 
the worldsheet curvature. The are the simple roots of the A r Lie algebra, (•, •) 
denotes the scalar product on the root space, p is the Weyl vector (half the sum 
of all positive roots) and the r-dimensional vector of fields can be expanded as 
= Y2i $i e %- The A r Toda theory is conformal provided Q and b are related via 



Q 



(2.2) 



The central charge is [20 



c = r + 12 Q 2 (p, p) = r 1 + (r + 1) (r + 2) [b + 



(2.3) 



The general form of a three-point correlation function in a 2d conformal field 
theory is [21] 



(V ai Zi) V Q2 (2:2, Z 2 )V a3 (z 3 , 3s)) 



C(ai,a 2 ,a 3 ) 



| /2l2 |2(A 1 +A 2 -A 3 )| 2l3 |2(A 1 +A3-A 2 )|^ 23 |2(A 2 +A3-A 1 ) ' 

(2.4) 

The Liouville theory is identical to the Ai Toda field theory and has a set of 
primary fields 



(2.5) 



The correlation function of three primary fields in the Liouville theory is [22] (see 
also K3I) 



C(a x , a 2 , a 3 ) — np>y(b )b 



-a 1 -a 2 -a 3 



X 



(2.6) 



x 



T(6)T(2 Ql )T(2a 2 )T(2a 3 ) 



T(«i + a 2 + « 3 — Q)T(— «i + a 2 + a 3 )T(ai — a 2 + «3)T(ai + a 2 — a 3 ) 



with 7 (6 2 ) = r ^ and T(z) = l/[T 2 (x|6, b~ l )Y 2 (Q - x\b, Zr 1 )] where Y 2 (z\e u e 2 ) 



is the Barnes double Gamma function 

In the Toda theories with r > 1 (W) primary fields can be defined in analogy 
with the Liouville case via 



(2.7) 



Recently it was shown [25, 26J that in the special case when one of the en's takes 
one of the two special values 



x = 



or 



x = *K , 



(2- 
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where Ai (A r ) is the highest weight of the fundamental (anti-fundamental) repre- 
sentation of the A r Lie algebra and x is a complex number, the three-point function 
is given by (12 .4p with 

(2Qp-a 1 -a 2 -X,p) 

C( ai ,a 2 , X ) = [n^(b 2 )b 2 - 2b2 \ ' x (2.9) 

x (W T(x) n e>0 T((Qp - a u e))T((Qp - a 2 , e)) 

Uij + («i - QP, h i) + ("2 - Qp, hj)) 

where the product in the numerator is over all positive roots and in the denominator 
the hi are the weights of the representation with highest weight Ai, cf. flA.2j) . (The 
result for x — x A r is obtained by replacing hi by h\ = —h r+2 -i-) 

Higher-point correlation functions in any CFT can be related to the three-point 
function of primary fields, which therefore determines the entire theory [2T]. Note 
that when r > 1 knowledge of the three-point function of W primary fields ( I2.9p 
does not determine all higher-point correlation functions, see e.g. [27| 15]. 

As an example, consider a four-point function. It is convenient to fix three points 
to 0, 1, oo and use a bra-ket notation which has the property (a\a) = 1, and is such 
that 

( ai |F a2 (l)y a3 (z)|a 4 > = (V Q - ai {0)V a2 {l)V a3 {z)V a4 {oo)) . (2.10) 
Inserting a complete set of states we find 

(ai \V a2 {l)V a3 {z) a 4 ) = / da > — — — — . (2.11) 

Here the intermediate states ^(c)) are descendants of the primary state labelled 
by a. (Throughout this paper we will label the internal momenta by a reserving the 
symbol a for the external momenta.) 

In the Liouville case it can be shown that (ipk>(a)\V a3 (z)\a4) is proportional to 
(cr|V^ 3 (;z)|a4) [21] and hence (12. lip can be calculated perturbatively. The ratio 

Ek^^il^li)!^^))^^)!^^^)' 1 ^^)!^^)!^) 

{ai\V a2 (l)\a){a\V a3 (z)\a A ) 

is called a conformal block. General n-point functions can be dealt with in an 
analogous manner. They depend on (n — 3) cross ratios. 

In the r > 1 case the situation is a little more involved, see [5] for a discussion. 



2.2 The A r quiver gauge theories and Nekrasov partition functions 

In jl] a class of conformal Ad Af — 2 generalised A r quiver gauge theories were 
introduced. This class of theories was denoted T( n>g )(A r ). The simplest example 
in this class of theories is the theory with a single SU(r + 1) gauge factor with 
2(r + 1) matter hypermultiplets in the fundamental representation of the gauge 
group. But the T( n , g )(A r ) class of theories includes many more theories, not all of 
which are conventional weakly-coupled gauge theories. The 7( n , ff )(A-) theories can 
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be viewed as arising from the six-dimensional A r (2,0) theory [2B] compactified on 
Cxi 4 where C is a genus g Riemann surface with n punctures. The genus of 
the Riemann surface depends on the number of loops in the (generalised) quiver 
diagram. The punctures are due to codimension 2 defects filling M 4 and intersecting 
C at points, and were argued in [3] to be classified by partitions of r + 1 (which 
can be represented graphically in terms of Young tableaux). One can therefore 
associate a Young tableau to each puncture. In the case of the A\ theories there 
is only one kind of non-trivial puncture. In the above example (a 7l,o{A r ) theory) 
there are two kinds of punctures. These are associated with the factors in the 
U(l) 2 SU(r + 1) 2 subgroup of the flavour symmetry group. The punctures associated 
with the SU(r + 1) factors are called full punctures and involve r mass parameters 
each and the punctures associated with the U(l) factors are called basic punctures 
and involve one mass parameter each. We refer to jl] for further details. 

A fundamental object in an M = 2 gauge theory is the Nekrasov partition 
function (from which the prepotential can be obtained). The partition function 
factorises into two parts as 

Z = Z pcrt Z inst , (2-13) 

where Z pert is the contribution from perturbative calculations (because of super- 
symmetry there are contributions only at tree and one-loop level), and Z- inst is the 
contribution from instantons. The most efficient method to obtain Z- inst is via the 
instanton counting method of Nekrasov [3]. This approach involves deforming the 
Af = 2 gauge theory with two parameters e\ and e 2 which belong to an SO (2) x SO (2) 
subgroup of the SO (4) Lorentz symmetry. We should stress that one needs the the- 
ory to be weakly coupled to be able to apply the instanton counting method. 

As an example, the instanton partition function in the SU(r + 1) theory with 
2(r + 1) fundamentals can be written [3] (see also [29] ) 

^-t = E^' ff U F( — v y v < 2 ' 14 ) 

_ , *- J \P>m Q>ni * mi * ni S)\E J \®"m Q>m *mj w 

y m,n=l sGY m v ' v v ' ' 

where y = e 2mT and the sum is over the (r+l)-dimensional vector of Young tableaux, 
Y = (Yi, Y2, . . . , Y r+ i) and \Y\ (the instanton number) is the total number of boxes 
in all the Y m 's. The a m parameterise the Coulomb branch of the theory and satisfy 
Yli=x ®i = ®- ft ^ S convenient to write a = J2i=i a i e « where ti are the simple roots of 
the A r Lie algebra. In the particular case of SU(2) this translates into a = (a, —a). 
In (T27T4|) e = ei + e 2 and 

E(x, Y m7 Y n , s)=x- e^yM + e 2 (A Y Js) + 1) , (2.15) 

where s = and i refers to the vertical position and j to the horizontal position 
of the box. Furthermore, Ly„ = fc n ,i — j an d Ay m = k^j — i, where k n> i is the length 
of the zth row of Y n and ■ is the height of the jth column of Y m . Finally, 

2r+2 

P(x, Yi, 8)=J[(x- (j - l)ei - {1 - l)e 2 - m f ) , (2.16) 

f=i 
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where the rrif are the masses of the matter fields (suitably defined). 

The perturbative (one-loop) piece in (I2.13p . Z pert , is a product of various factors. 
For SU(r + 1) the gauge field contributes a factor 



II Fr -~ 1 vr~r « 1 V ( 2 - 17 ) 

- LJ - 1 2 {ai - % - e 2 |ei,e 2 JI 2 (aj - % - ei|ei,e 2 J 



where r 2 (a;|ei, e 2 ) is the Barnes double gamma function [21], and each of the massive 
hypermultiplets transforming in the fundamental representation of the gauge group 
contributes a factor 

r 

J]r 2 (a i -m / + e|e 1 ,e 2 ). (2.18) 



2.3 The AGT relation 

The AGT relation is a relation between the two classes of theories discussed in 
sections 12.11 and 12 .21 Up to an overall factor, the instanton partition function of 
a T(n, g )(A r ) theory is (conjectured to be) equal to a chiral block of an n-point 
correlation function in the A r Toda field theory formulated on a genus g surface. In 
other words, the punctures correspond to insertions of vertex operators in the Toda 
theory. For example, the instanton partition function in the SU(r + 1) theory with 
2r + 2 fundamentals is equal to the chiral block (in a specific channel) of the four- 
point function in the A r Toda theory on the sphere. The momenta of the vertex 
operators, a i: are mapped to the masses, m ; , in the gauge theory. This relation 
is linear (the exact form depends on conventions for the gauge theory masses). 
Furthermore, the internal momenta in the chiral block, <7fc, are linearly related to 
the dk Coulomb moduli. Finally, the parameters e\ and e 2 in the instanton partition 
function are related to the parameter b in the Toda theory via, 

The most common choice is to set e\ — b, e 2 = |. There is also a slight extension of 
the above result where the full partition function including the perturbative piece is 
related to the full correlation function including the three-point pieces. For further 
details about the AGT relation, see [2j [5] . 



3 The matrix model approach of Dijkgraaf and Vafa 

The first step in the analysis of [1] is to realise the relevant quiver gauge theories 
in string theory using geometric engineering [15] and then use the results in [30] 
to relate the topological string partition function (which is equal to the Nekrasov 
partition function) to a matrix model calculation. This argument works provided 
ei = — e 2 = g s . How to deal with the case of general e ij2 was also proposed in [TJ 
and will be discussed later in this section. 
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In the case of interest to us the relevant geometries are A r singularities 



uv — X 



r+1 



0. 



(3-1) 



(There is also a decoupled C factor parametrised by z.) The matrix model corre- 
sponding to this geometry is the so called A r quiver matrix model [THUS]- It involves 
r hermitian iVj x iVj matrices $j (i — 1, 
The matrix model partition function is (proportional to) 



, r) as well as the N{ x Nj matrices Bij. 



exp 



-- JJtr [(2% - AiMBi&Bji - Bj&iBit)] 



i<j 



(3.2) 



where A^- is the Cartan matrix of the A r Lie algebra, i.e., = (e$, e^). In terms of 
the eigenvalues of $j, \{ (1 = 1,..., iVj), the partition function (13. 2 j) becomes, after 
integrating out the Sjj, 



]JdXl J] |Af - A 

(i,I)<U,J) 



J\A Z 



(3.3) 



Here (i, I) < (j, J) if i < j or i = j, I < J. 

There is a close relationship between the A r matrix models and conformal field 
theory [T6l \T7\. In this relation r free bosons, (fi, in a 2d conformal field theory are 
related to matrix model quantities in the following way^ 



dtpi(z) = — tr 
9s 



(3.4) 



Because of the relation (13.41) the free boson CFT vertex operator V a (z) = e a '^'^ 
translates into 

JJdet(2-$i) Ql/9s . (3.5) 

i 

Here ot % = (ei,a) where a = ^a'Aj (see appendix IA.1I for a summary of the Lie 
algebra terminology). Correlation functions of (free boson) CFT vertex operators 
can therefore be calculated using matrix model technology. The correlation functional 



(V ai /g s {zi) ■ ■ ■ V ak / gs (zk)) ao ,N 



(3.6) 



translates into 



/ n d$ * n dB ^ n det ^ - ^ ai/gs n det ^ - ^ 

— JJtr ((Sij - Aij)(B ij § J B ji - W'Vti 



>ai/g„ 



(3.7) 



2 This relation holds when the matrix model potential is zero which is the case we are interested 
in. The factor of g s is unconventional but convenient for our purposes. 

3 In this paper we restrict our attention to correlation functions on genus surfaces. 
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In (13. 6p the subscript «o refers to the fact that there is an extra olq charge at infinity. 
This means that the correlation function really involves k + 1 vertex operators. 
Furthermore, in our conventions, 

k 

A l] N ] =a j Q /g s -Y,<l9s- (3.8) 

n=l 

To analyze the matrix model expression (13 .7\\ one can use the following relation 
det(z — $) a / 9s = exp(^-tr \og(z — $) ). The insertions therefore effectively induce 
the matrix model potentials (of multi Penner type) 

k 

Wii^i) = tr ^o£]og(2 a - $0 , (3.9) 

a=l 

and the matrix models can therefore be analyzed using standard techniques. 
In terms of the eigenvalues the matrix model correlation function is 

/ ri dA ' n i a ' - a /i Aij - xi ^ ai/9s ■ ■ ■ n& - x ^ ai/9s - ^ 

il {i,I)<{j,J) i,I i,J 

In this expression we have left the integration contour unspecified. The choice 
of integration contour turns out to be quite subtle and will be discussed in later 
sections. In previous applications of matrix models to supersymmetric gauge theories 
[30J the matrix models should properly be thought of as holomorphic matrix models 
with a choice of contour, see, e.g., [32] for a discussion. In the present case, there 
are also additional subtleties since the potentials are logarithmic. 

The proposal of Dijkgraaf and Vafa [1] is as follows: to connect the chiral corre- 
lation functions of the vertex operators V involving the free scalar fields (fii to the 
correlation functions of the chiral vertex operators V involving the fields 0, in the A r 
Toda field theory one should take the large N limit and identify the a a 'a (including 
«o) with the external momenta of the Toda theory vertex operators. Furthermore, 
the matrix model potential (I3.9P has stationary points which in the large N limit ex- 
pand into cuts. The corresponding filling fractions gsNf 1 (subject to the constraint 
Yl^=i = Ni) are related to the internal momenta a m in the Toda theory (m 
label the internal momenta and % label the components of each of the u m ). 

This proposal shares many similarities with the earlier work |30j, but we should 
stress that in [I] the number of terms in the potential is related to the number of 
nodes in the gauge theory quiver, whereas the number of matrices is related to the 
rank of the gauge group; in |30j the roles were reversed. 

In [Tj there is also a discussion of the AGT relation using brane probes; this 
approach will not be used in this paper. 

The analysis so far only involves g s , i.e., e\ = —62. As mentioned above there 
is a further refinement of the matrix model that is needed to treat the case with 
general ei^. In [TJ it was suggested that the required modification is the so called (3 
deformation (or /3 ensemble) [52] ■ This deformation changes (I3.10p to 

/ri dA ' n iA^-A/i^n^- A ')^ Qi/9s ---n^- A ')^ /9s ( 3 - n ) 

il (i,I)<(j,J) i,I i,I 

8 



with the identification j3 = — e 2 /ei and g s = \/— eie 2 . Also, ( 13. 8 p changes to 

0A .. Nj = ^l a ^^lj2<. (3.12) 

9s 9s n=l 

We should stress that for general j3 the above integral (13. lip can no longer be viewed 
as arising from an integral over matrices in any reasonable way. Therefore, strictly 
speaking, we are no longer dealing with a matrix model. Sometimes the model for 
general (3 is called a generalised matrix model, but we will by a slight abuse of 
terminology continue to call it a matrix model. 

In the next section we analyze various aspects of the above matrix model for 
the case of the A\ theory and make detailed calculations and comparisons with the 
corresponding expressions in the 4d Ai quiver gauge theories and the 2d Liouville 
theory. In section a similar analysis will be performed for the A r theories with 
r > 1. 



4 The A\ matrix model 

In this section we perform several calculations in the A\ matrix model. The resulting 
expressions are compared to the corresponding expressions in the Liouville theory 
and the A 1 quiver gauge theories. 



4.1 The three-point function 

Our first example is the matrix model three-point function: 



KJ I 



Note that l/e x = y/]5/gH This complicated looking integral is the so-called Selberg 
integral [M] which can be evaluated exactly with the resultEl 



N 



1 r(2a 1 /e 1 + ! + (/- l)f3)Y{2a 2 /e 1 + ! + (/- l)P)T(ip) 

(27r)^Al T(2a 1 /s 1 + 2a 2 /s 1 + 2 + (I + N-2)^T((3) ' [ ' } 

In evaluating the above integral we assumed that the choice of integration contour 
is such that the A 7 's are integrated over the interval [0, 1]. (It is possible to perform 
changes of variables in the above integral to obtain other integration ranges; see 
e.g. [33j section 17.5].) Using the result 

m = v^(- g ,) 1/2 - r ^ ( "" ll 7"" £2) ; . ("-si 

1 2\-ze\ - e 2 \ - ei, -e 2 ) 



4 When referring to matrix model quantities we will use the notation rather than since it 
will turn out that our conventions are such that Ej = — e*. 

5 In appendix lA.2l we present an alternative derivation of this result when /3 = 1 using orthogonal 
polynomials. 
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where T 2 (x| — £1,-62) is the Barnes double gamma function [24J, together with 



P = - a and 



Nz T 2 (-ze x + {N-l)e 2 \ 


- £l, -£2) 


r 2 (-z£i - £ 2 


-£l,-£ 2 ) 



7=1 

we find that ( 14. 2 p equals (here T 2 (x) is short-hand for r 2 (x| — £1, —£2) and e = £i+£ 2 ) 

%r) x (4 - 5) 

r 2 (-2ai + Ne 2 - e)T 2 (-2a 2 + Ne 2 - e)T 2 (N e 2 )T 2 (-2a 1 - 2a 2 + Ne 2 - 2e) 
r 2 (-2ai - £)r 2 (-2a 2 - £)r 2 (0)r 2 (-2ai - 2« 2 + 2iV£ 2 - 2e) 

Finally, using Ne 2 = —N(3ei = (— cto + a.\ + a 2 ) we obtain 

r 2 (-ao-«i+«2-£)r 2 (-a +ai-a 2 -£)r2(-ao+ai+«2)r2(-a -ai-Q!2-2£) 
r 2 (-2ai - £)r 2 (-2a 2 - £)r 2 (0)r 2 (-2a - 2e) ' 

In general, the three-point function in a 2d conformal field theory does not factorise 
into holomorphic and anti-holomorphic parts so there is no unambiguous meaning to 
a 'chiral three-point function'. However, after suitably rescaling the vertex operators 
with multiplicative factors depending on their momenta, the Liouville three-point 
function (12. 6p can be written as (recall that a* = Q — a^) 

[T(ai + a 2 + a 3 - Q)T(— «i + a 2 + a 3 )T(ai - a 2 + a 3 )T(ai + a 2 - a 3 )] _1 = 
\T b (2Q-ai-a 2 ~a 3 )T b (Q+a 1 -a 2 -as)T 2 (Q-a 1 +a 2 -a 3 )T 2 (a 1 +a 2 -a 3 )\ 2 , (4.7) 

where r&(x) is short-hand for T 2 (x\b, b^ 1 ). Therefore, for the Liouville theory there 
is a natural definition of a chiral three-point function as the "square root" of (14.71) . 
After a suitable redefinition of the matrix model vertex operators, we see that the 
matrix model three-point function (14.61) precisely captures the chiral part of the 
Liouville three-point function, provided that «o —> a 3 and we identify 

ei = -b, £ 2 = -1/6. (4.8) 

The matrix model expression can also be compared with (the perturbative part 
of) the Nekrasov partition function of the corresponding gauge theory, which in 
the present case is the so called T 2 (or 7^,0(^1)) theory — a theory of four free 
hypermultiplets |4j. Redefining the matrix model vertex operators as above we are 
left with the four T 2 factors in the numerator of (14.61) . These are of precisely the 
right form to reproduce (the perturbative part of) the T 2 theory, cf. (12.181) (with a 
suitable definition of the four masses). (Possibly one can also make sense of ( 14. 6 p 
within the framework in [35].) 

We should also mention that the derivation in |22j of the full three-point function 
in the Liouville theory was based on an argument which involved a complex version 
of the above integral. Although the equations are similar the logic was different. 
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4.2 Higher— point functions 

One can also consider higher-point functions in the matrix model. A tractable 
example is a four-point function, where one of the «j, a 3 say, is equal to E\j2 or 
£2/2, i.e. the integral 



/ ri dA/ n i y - xJ \ w - a 1 ) 2 ^^) 2 ^! - ^f a2/£i 
1 i<j 1 



(4.9) 



with «3 equal to e±/2 or £2/2. In [36J it was shown that the above integral satisfies 
the hypergeometric differential equation: 

z{l-z) ^fi- + [C-(A + B + l)z] — AB F{z) = , (4.10) 
where, if a 3 = £i/2, 

A = -N, B=-(2a 1 /ei + 2a 2 /ei + 2) + N-l, C = ^{2a 1 /e 1 + 1) , (4.11) 
and, if 03 = £2/2, 

A = (3N, B = -(2a 1 /ei+2a 2 /e 1 + l) + (3(2- N) , C = -2a x /ei + /3 . (4.12) 

On the other hand, it is known that in the Liouville theory the four-point func- 
tion 

(V. b/2 (z)V ai (0)V a ,(l)V ao (oo)) , (4.13) 

satisfies the equation 



-b 



z{z-l) * z 2 {z-V, 

2 A(-6/2) + A( ai ) + A(a 2 ) - A(a ) 



iZ"(z) = 0, (4.14) 



where A (a) = a(Q — a)- After writing H(z) = z bai (l — z) boi2 F(z) the above equation 
reduces to the hypergeometric equation (14. 10p with 

3 1 

A = b(a 1 + a 2 -a -b/2) , B = b(a x + a 2 + a - -b- -) , C = 6(2«i -b) . (4.15) 

2 

Using /3 = — ^ and —e 2 N — a + a.\ + «2 + £i/2 = (compared to the corresponding 
expression for the three-point function there is now an extra term coming from 
«3 = £i/2), we see that (14. lip and (I4.15P agree provided we use the identifications 
(14.81) . This analysis shows that the matrix model integral (14. 9 j) with the above 
identifications is proportional to the chiral block in the Liouville CFT. Note that for 
the special correlation function (14. 13ft the internal momentum is restricted to two 
discrete values corresponding to the two solutions to the hypergeometric equation. 
The case with an insertion of the vertex operator V_i(z) can be treated analogously 

b 
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and compared to ( I4.12p . The choice of signs in ( 14. 8 p differs from the usual convention 
and means that = — but note that both eie 2 and ^ are unaffected. We could 
make sign changes elsewhere to restore the usual rule but this would clutter some 
of the above formulae. 

One can also check that the above expression agrees with the corresponding 
Nekrasov partition function. This was anticipated in [5] and derived in detail in [6] . 
We briefly recall the argument here. The instanton partition function is as in (J47T7J) 
with m, n — 1, 2 and (Si, d 2 ) = {a, —a). If we tune the Coulomb modulus a to fulfill 
P{a) = by setting a = mi, then only terms in the sum with only Y 2 non-empty 
give non-vanishing contributions. If we furthermore set m 2 = —a — e\ then only 
those Y 2 tableaux that have boxes only in the first column survive (in other words 
k 2 j is only non-zero for j — 1, so that i — 1, . . . , k 2l and k 2 ^ = 1)- Next using the 
AGT relation in the form 

mi = -^+tti+a 3 , m 2 = ^-«i+a 3 , m 3 = -|+a 4 +a 2 , rn 4 = ^-a 4 +a 2 , (4.16) 
one finds a 3 = — ei/2 and 

^ (A)i(BW Ul7) 

^t-^-^^, (4.17) 

where (X) n = X(X + 1) • • • (X + n — 1) is the Pochhammer symbol and 

3 

A = (ai+a 2 -a 4 -y)/e 2 , B = {ai+a 2 +a i --e 1 -e 2 ) /e 2 , C = -(2a 1 -e 1 )/e 2 , 

(4.18) 

which agrees with (14. 1 5j) provided a 4 — > a , e\ = b and e 2 = \. Together with the 
fact that the expression (14.171) is precisely the series expansion of the hypergeometric 
function 2 Fi(A, B; C;y), which solves the differential equation (I4.10p with y = z, 
this shows that the instanton partition function agrees with the chiral block (up to 
an overall factor). 

The matrix integral corresponding to a correlation function with k insertions of 
V-bfi i-e- 

k 

J I KJ I a=l 

was also calculated exactly in |36j. The result is (proportional to) the generalised 
hypergeometric functional 

2 if (-N, h[2a x + 2a 2 )/e 1 + k + 1) + N - 1; ^{2a 1 /e 1 + k); z u . . . , z k ) . (4.20) 
The function 2 F^(Ai, A 2 , B\ \ Zi, . . . , z k ) is defined as 

2 F?(A, B; C; Zl ,..., z k ) = £ ^ ^ 1 , (4.21) 



6 Note that it is also possible treat the cases with V-u,l2 insertions by setting z\ 
dUgi and (TODl) . 
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where the sum is over all partitions £ with at most k parts, 

= IK* - ^ - > ( 4 - 22 ) 

and Pf(zi, . . . , Zk) is a (properly normalised) Jack polynomial. See [36] for further 
details about the notation. It should be possible to show that the corresponding 
CFT and gauge theory calculations lead to the same result. 

We should also point out that the treatment of gauge theory surface and line 
operators in the Liouville language discussed in [9j involves the insertion of vertex 
operators with a = —6/2. This is precisely the class of operator insertions we have 
discussed above. 

To go beyond the restricted set of correlation functions discussed above one 
possible approach would be to try to mimic the CFT method and determine the 
correlation functions in a perturbative expansion in the £j. 

An alternative approach is to use matrix model perturbation theory. For the 
four-point function one can (at least when (3 = 1) use the method in [371 EH]; f° r the 
multi-cut and A r quiver extensions, see e.g. [3U]. One can also obtain a perturbative 
expansion within the framework of sect ions [4 . 3 1 and l4~4l see section l4"751 below for some 
sample calculations. On the CFT side the matrix model perturbation theory is a 
somewhat peculiar expansion and it is not clear what its relevance is. A drawback of 
the perturbative matrix model approach is that solving for the stationary points lead 
to complicated expressions, but a clear advantage is that one can handle arbitrary 
punctures {i.e. arbitrary vertex operator momenta) using this method (both for the 
A\ and A r theories). Let us also mention that in [2H] another, less direct, matrix 
model approach to Af = 2 gauge theories was discussed; it might be interesting to 
try to connect it to the present approach. 

4.3 The curve: one— cut solutions 

In this section we analyse the one cut matrix model spectral curves, focusing on the 
matrix model corresponding to the three-point function discussed in section |4~T1 We 
start by reviewing some well-known results. See e.g. [ID] for further details. In the 
standard diagonal gauge (eigenvalue basis) the one-matrix model partition function 

Z = * [f]dX I T](X I -X J ) 2 exp(--yW(X 1 )] 

(2nf N\J l\ }}; [ { 9s jrt 7 

1 f N 

= {2n)N m J U dXl ^V(N 2 S CS [p(A)]), (4.23) 

where 

S eS \p(X)] = -j [ dAp(A)^(A)+ // dAdA / p(A)p(A / )log|A-A / | . (4.24) 

t JC J JCxC 

To conform with standard matrix model conventions W in this section and in sections 14.41 and 
14.51 corresponds to — W elsewhere in the paper. This also implies that the ctj differ by a sign. 
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Here t is the 't Hooft coupling t = Ng s and we have introduced the eigenvalue 
density p(X) = J2f=i ~~ ^ £ )i normalized as J c dX p(A) = 1. In this expression 
one still needs to specify the geometrical nature of the cut, C. In the most general 
case p(A) has compact support, with C a multi-cut region with s cuts. For the 
moment we shall focus on the one cut case, with C = [a,b]. If one now considers 
the Riemann surface which corresponds to a double-sheet covering of the complex 
plane, C, with precisely the above cut, it is natural to define the A-cycle as the 
cycle around the cut. In this case, the -B-cycle goes from the endpoint of the cut to 
infinity on one of the two sheets and back again on the other. 

The generator of single-trace correlation functions is given by the resolvent 

i / i \ i +°° i 

" w= if ^b^ 1 ^' (4 - 25) 

x ' fc=0 

which has the standard expansion u(z) = Ylg^o 9"s 9uj g( z ) with 

u (z) = J dX-^-. (4.26) 

The normalization of the eigenvalue density then implies that uq{z) ~ ~ as z — > +oo. 
Also, observe that ojq{z) is singular for z G C while it is analytic for z C. One 
may compute coo(z) by making use of the large N saddle-point equations of motion 
of the matrix model, 

u (z + ie) + u (z - ie) = -W'(z) = 2 PV f dX (4.27) 

t Jc z- X 

In a similar fashion, coq(z) is related to the eigenvalue density as 

p(z) = -— (uj (z + ie) - u {z - ie)) = — hi^(z). (4.28) 

Z7T1 7T 

For a generic one-cut solution, the large resolvent is given by the ansatz 



^.l^gjt'-f-') (4.29) 
2t J c 2ni z — w V (w — a){w — b) 

where one still needs to specify the endpoints of the cut, {a, b}. An equivalent way 
to describe the matrix model geometry is via the corresponding spectral curve, y(z), 
which basically describes the geometry of the Riemann surface we mentioned above. 
One may write 

y ( z ) = W\z) - 2t u (z) = M(z) y/ ( z - a)(z - b), (4.30) 

wittH 

= r a. iniM i (4 . 31) 

/ (0) 27Ti 1-wz y (1 — aw)(l — bw) 



3 This particular expression only holds for polynomial potentials. 
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where, again, one needs to specify the endpoints of the cut, {a, b}. The aforemen- 
tioned large z asymptotics of the resolvent immediately yield 2 conditions for these 
2 unknowns. They are 

dW W " W ' M 2U m (4.32) 



2vri y/( w - a )(w - b) 

for n — 0, 1, fully determining the endpoints of the single cut s — 1. 

It is also useful to define the holomorphic effective potential as V^. eS {z) = y(z). 
The effective potential is then given by the real part of the holomorphic effective 
potential, in such a way that 

Kff(A) =Re j dzy(z). (4.33) 

J a 

The real part of the spectral curve therefore corresponds to the force exerted on a 
given eigenvalue. The imaginary part of the spectral curve, on the other hand, is 
related to the eigenvalue density via 

p{z) = ^-lmy{z). (4.34) 

Finally, it turns out that one may also write the 't Hooft parameter in terms of the 
spectral geometry as 

t = -^j> dzy{z). (4.35) 

We now turn to our main point and consider the large N expansion of the matrix 
model with potential 

k 

W(z)=^2on\og(z-Zi) (4.36) 
i=i 

where k, {a:j}* =1 and {zj}f =1 are parameters we shall keep unspecified for the mo- 
ment. It follows from ( I4.36P that 

1 = 1 1 

implying that the logarithmic terms in f 14 . 3 6 1) will not be terribly problematic — one 
only needs to take into account extra poles, when moving the contours of integration 
around the complex plane. We begin by focusing on the one-cut solution, for which 
the large N resolvent is given by the ansatz 



1 fdwW'iw) (z-a)(z-b) , AOO , 

u Q (z) = — <b —\t 77 rV, (4-38) 

v; 2t J c 2vri z-w V (w-a)(w-b) V ; 

where the integrand now has poles at the locations {zi} but, because the potential 
is purely logarithmic, any pole of the integrand at infinity is gone. In this 
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straightforward deformation of the integration contour reduces the integral along 
the cut to a sum of simple poles as 



i (z - Zi) yj(zi - a){zi - b) 



y/(z-a)(z-b) . (4.39) 



The large z asymptotics, co (z) ~ - + • • • as z — > oo, immediately implies that the 
endpoints of the cut C = [a, b] are determined by the system 



E 

i- 



2<%- 



■( yj(zi - a){zi - b) 

lOLiZi 



2t, 



y/(zi - a){zi - b) 

and the single-cut spectral geometry is then described by the curve 

2a, 



(4.40) 
(4.41) 



E 



i (z- z^ y/(zi- a)(zi-b) 



y/(z- a)(z - b). 



(4.42) 



One may also compute the holomorphic effective potential in a simple manner. We 
obtain 



V h - eS (z) = -2 ^2t-J2 2a * ) log 2 (v/F^ + v^fc 
+ ^ 2^10, 



8=1 



y/z — a y/z~i — b 



y/z -by/% 



2a.i log 



+ 



1 + 



(4.43) 



yj Z — ay/Zi — b 

y/z — b y/Zi — a 



with a and b determined by the system above. The structure of Stokes lines for this 
effective potential will be more complicated than in the usual polynomial cases. 

Let us now specialise to the matrix model corresponding to the three-point 
function, with potential 



Wapf (z) = 2cui log z + 2a 2 log (z - 1) , 
and the further constraint 



t = a + 



0£i. 



(4.44) 



(4.45) 



i=i 



Next we turn to the study of the spectral geometry associated with the matrix 
model potential ( I4.44p beginning with the "classical" geometry. The critical points 
are located at the points z* such that 



. , . 2ct\ 2ao 
W^(z^) = ^ + L- = 0. 



(4.46) 
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In this case, the general solution is 



z* = (4.47) 
«i + a 2 

In the classical limit where the 't Hooft coupling vanishes (i.e., where we choose 
vertex operators such that ao + ol\ + a 2 = 0), one simply has y(z) = W^ p{ (z) and 
the cut collapses to the critical point of the potential, [a, b] — >■ z m . In this particular 

case, 

Wspf = 2a, log — ^— + 2a 2 log ( 4 .48) 

Now, the full spectral geometry will be such that the critical point z* opens up into 
a branch cut, of size t (and not touching the marked points associated with the 
vertex operator insertions at {0,1, oo}). This blown-up geometry of the spectral 
curve will have its shape determined by the parameters «o, a, and a 2 . Clearly, 
because there is a single critical point, the spectral geometry will correspondingly 
have a single cut — the situation we studied above. The spectral geometry associated 
to the matrix model with potential (I4.44p is therefore a genus zero one-cut Riemann 
surface. From our previous general results it follows that the spectral curve is 

y(z) = [ -^L + 2a2 ) y/( z - a )(z-b), (4.49) 

V; \z^b (^-1)^/(1-^(1-6) J VV V ^ 

where the endpoints a and b are obtained from the solution to the system 



2a, , 2a 2 = 



ab - a)(l - b) 

2a 2 

y/(l-a)(l-b) 



-2a , (4.50) 



with (partial) solution 



ab = — , (4.51) 
a 
a 2 



y/(l-a)(l-b) = - — . (4.52) 

a 

This immediately simplifies the spectral curve to 

V(z) = - ( ^_° z) y/(z-a)(z-b), (4.53) 
For completeness we also give the explicit solution to f)4.50p : 

al+af— a\— a/(« — a,— a 2 )(a +ai— a 2 )(a Q — ai+a 2 )(a +ai+a 2 ) 



2a\ 



( al+al—al+^(a —ai—a 2 )(a +ai—a 2 )(a —ai+a 2 )(a +ai+a 2 ] 
b = A L. (4.54) 



2ai 
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The spectral curve ( I4.53P can be written 

z 2 (l-z 



v 2 = r£^, (4.55) 



where P2(z) is a polynomial of degree 2. This expression was also written in [T]; 
here we have also explicitly determined the coefficients in Pi{z) in terms of the three 

oti's. 



4.4 The curve: multi-cut solutions 

Let us return to the geometrical nature of the cut, C, where C is now a multi-cut 
region with s cuts. There are two cases: when s is smaller or equal to the number 
of minima of the potential V(A), a typical situation in the standard matrix model 
context; or when s is equal to the number of non-degenerate extrema of the potential 
V(A), the situation which arises when dealing with topological strings which is the 
case relevant to us. More precisely 



C = |J A 1 , (4.56) 



i=i 



where A 1 = [X21-1, are the s cuts and x\ < X2 < ■ ■ ■ < X2 S - If one now considers 
the hyperelliptic Riemann surface which corresponds to a double-sheet covering of 
the complex plane, C, with precisely the same cuts as above, A 1 , it is then natural 
to define the v4 7 -cycle as the cycle around the A 1 cut, with the i?/- cycle following 
via Bj fl A J = 5j. In this case, the S/-cycle goes from the endpoint of the A 1 cut 
to infinity on one of the two sheets and back again on the other. 

For a generic multi-cut solution, the large N resolvent is given by the ansatz 



1 f dw W'(w 



2s 

z-x k 



11^-^ (4.57) 



2t J c 2m z — w \^ jJ-Jj w — Xk 

where one still needs to specify the endpoints of the s cuts, {xk}- An equivalent 
way to describe the matrix model geometry is via the corresponding spectral curve, 
y(z), which basically describes the hyperelliptic geometry of the Riemann surface 
we mentioned above. One may write 

y ( z ) = W\z) - 2too (z) = M(z) v/^sM, (4.58) 

where 

2s 

Y[(z-x k ) (4.59) 



2s 
k=l 



and 9 ' 



M(z) = <b y - 1 -^- , (4.60) 

l(n\ 2m 1 — WZ /tt2s / 1 \ 

(0) \ Uk^i 1 - x kW) 



'This particular expression only holds for polynomial potentials. 
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and where, again, one still needs to specify the endpoints of the s cuts, {x^}. The 
aforementioned large z asymptotics of the resolvent immediately yield s + 1 condi- 
tions for these 2s unknowns. They are 

dw w n W'(w) n . . 

— = 2t5 ns , (4.61) 



2m \/U 2 k=i(w-x k ) 

for n = 0,1,. . . , s. In order to fully solve the problem, one still requires s — 1 extra 
conditional for the full set of {xk}- These extra conditions depend on whether 
one wants to consider the standard matrix model or the topological string case. In 
the first option one considers all the different cuts at equipotential lines, where this 
condition may be written as 

dzy(z) = 0. (4.62) 

A physical understanding of this expression says that there is no force moving eigen- 
values from one cut to another. In contrast, the topological string option (which is 
the case relevant to our analysis) generically corresponds to an unstable situation 
from a purely matrix model point of view. In this case one considers the filling 
fractions, 

at r 

V* = jf = J dX PW> I =1,2,..., s, (4.63) 

as parameters, or moduli, of the problem under consideration. Observe that here 
^j =1 T] 1 = 1, making it an actual total of s — 1 extra parameters, precisely the 
number required. By re- writing the eigenvalue density in terms of the resolvent, and 
the resolvent in terms of the spectral curve, one is led to the equivalent definition 



r] 1 = -^-r <f> dzy(z). (4.64) 



Airit 



One may also use as moduli the partial 't Hooft couplings t 1 = tr\ l = g s Nj. In this 
case 

t 1 = — <p dzy(z), (4.65) 



47ri 



A 1 



with X)/=i — t, making a total of s — 1 moduli. 

Let us now return to the large N expansion of the matrix model with the potential 
(I4.36P and associated derivative (I4.37P and briefly discuss how one may use standard 
saddle-point techniques to address multi-cut solutions. Again the logarithmic terms 
are not a problem; as we have seen before one only needs to take into account 
extra poles, when moving around the complex plane. When addressing multi-cut 
solutions, with s cuts, the large N resolvent is given by the ansatz (14.571) . where the 
integrand now has poles at the locations Zi and any pole of the integrand at infinity 



10 Observe that no further conditions were required in the previous one-cut case, where C = 
[xifXz] = [a,b]. Indeed, in that situation the large z asymptotics fully determined the endpoints 
of the single cut. 
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is gone. As in the single-cut straightforward deformation of the integration 

contour reduces the integral along the cut to a sum of simple poles as 



u (z) = i [ W\z) - 



2ai 



i=l [Z — Zi 



Wn=\ ( Z i X n) 



2s 



\ n=l 



(4.66) 



The large z asymptotics, ojq(z) 



as z — > oo, immediately yield s + 1 



conditions on the 2s endpoints of the cuts via the system 
k 



E 

i=l 



2a,- vf 



dw 

C 2m ( W ~ Z *) \lYFn=l( W - X n) 



2t 5 r , 



(4.67) 



for m = 0,1, ... ,s. This may be written more explicitly as (notice that now the 
integrand does have a pole at infinity, when m = s) 



E 



E 



2a, ; 



1 \/rU=i (.Zi-Xn) 



1 \/rL=l (Zi-Xn) 



E 



E \ 2a - 



i=i 



2a,- z* 



s-l 



2a» zf 

rin=l ( Z i ~ X n) , 



It. 



(4.68) 



Of course in order to fully solve the problem one still requires s — 1 extra conditions 
for the full set of endpoints {xk}- Finally, the spectral geometry is described by the 
hyperelliptic spectral curve 



E 



2ai 



i=l [Z — Zj 



2s 



\ n=l 



(4.69) 



Having understood the multi-cut spectral geometry we now focus on the case 
corresponding to the chiral four-point function in the Liouville theory, namely the 
matrix model with potential"] 



W ip{ (z) = 2«i log z + 2a 2 log (z - 1) + 2a 3 log (z - () . 



(4.70) 



11 Here and in section 03] for clarity we use ( = z\ to denote the location of the vertex operator 
insertion; in other sections z is used. 
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We begin with the "classical" geometry of the potential. The critical points are 
located at the points z* such that 

In this case, the general solutions are 



(1+C) «i+C «2+«3- V ((1+0 ai+C "2+«3) — 4ai C (ai+a 2 +a3) 

2 (ai + a 2 + a 3 ) 



(1+C) ai+C a 2 +«3+\/ ((1 + "1+C tt2+«3) 2 — 4ai C (ai+a2+«3) 

<z*, 2 = ^7 x -(4-72) 

2 (« 1 +a 2 +o;3j 

Generically, the critical points z^x and z*^ will open up into branch cuts, of sizes t\ 
and t% (and not touching the marked points associated to the vertex operator inser- 
tions at {0, 1, 00}). Because of the two critical points, the most general spectral 
geometry will correspondingly have two-cuts and the spectral geometry associated 
with the chiral Liouville four-point function is a genus one two-cut (elliptic) Rie- 
mann surface. This two-cut blown-up geometry of the spectral curve will have its 
shape determined by «o, «i, a 2 and a^. To be more precise, from the large z asymp- 
totics of the genus zero resolvent one obtains 3 conditions on the endpoints of the 
two cuts; the remaining required condition arising from the partial 't Hooft moduli, 
ti or t 2 (where t\ +t 2 = t). In [T] this modulus is actually traded for a = t 2 — ti, the 
Coulomb modulus in the gauge theory, and we shall use this notation henceforth. 

Of course there are particular points in the moduli space of the elliptic spectral 
curve where the geometry simplifies. One is the degenerate case where both partial 
't Hooft couplings vanish. Another special point occurs when only one of the critical 
points opens up into a branch cut, in which case one is dealing with a one-cut 
pinched spectral geometry, the pinch at the location of the critical point that remains 
"closed" . Let us consider this special case, where 



I a + i^ a . = , (4.73) 

i=0 

I 1 3 

a+ J2^ = t. (4.74) 



2 2 

i=0 



From our previous (single-cut) result (I4.42p we have 
y(z) = y/{z- a)(z - b) 



x 



2a\ 2a 2 2a^ 

+ ; + 



ab (z-l)y/(l-a)(l-b) (z-0 y/(C-a)(C-b) 



,(4.75) 
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where the endpoints a and b solution to the system 



2ai ' 2Q2 2a3 0, (4.76) 



ab y/(l-a)(l-b) y/(C-a)(C-b) 

+ s = -2a . (4.77) 



-a)(l-6) >/(C-a)(C-6) 
This may be equivalently written as 



2ai 2. 0+ / C - 1)2Q3 , (4.78) 



~ 2a o m w , =t. ( 4 - 79 ) 



which slightly simplifies the spectral curve to 

y(z) = f-^g- + s ; 1 y/ (,-a)(,-&) . (4.80) 

Generically, however, we are dealing with a system of quartic equations and, al- 
though it can be solved algebraically, its exact solution is not terribly illuminating. 
In the following, we therefore choose a different route and solve this system pertur- 
batively in (. This is motivated by the expansion on the CFT side and simplifies 
the problem considerably. To first order, we obtain the solution 

= ft i + »3 _ ((<*! + «3) 2 + ol\ - og) a 3 ( Q / 3/2 v 
"o 2a (ofi + « 3 ) 2 

V(l "«)(!-») = f 1 " + O (C 2 )) , (4.81) 

or, equivalently, 

a = 2% + 2ag(a 1 + a 3 ) <^ + (0. 

2ag 2a;5 («i + a 3 ) 

where we have defined 

fi = (ao+ai+«2+«3) («o— ai+a 2 — a 3 ) (a +ai— a 2 +a 3 ) (a — a\— a 2 — a 3 ) . (4.83) 

This solution will be important in the perturbative calculations in section l4"75l below. 

From having worked out some degenerate cases of the matrix model associated 
with the Liouville four-point function we have acquired some intuition about what 
to expect as we move on to the case with arbitrary vertex operators with conformal 
dimensions such that, generically, one will find a two-cut geometry. Let us briefly 



22 



comment on this geometry. We have previously computed the spectral curve for 
a general multi-cut case and, for the two-cut ansatz associated to the Liouville 
four-point function, this is 



2ot\ 2a<2 



z y /x 1 x 2 x 3 x 4 (z - 1) a/(1 - Xi) (1 - x 2 ) (1 - x 3 ) (1 - X A ) 
+ ) (4-84) 



with 



c(z) = (z - Xi) (z - x 2 ) (z - x 3 ) (z - x A ) . (4.85) 

It is simple to see that the endpoints of the two cuts, {xi}f =1 , are now a solution to 
the system 

= (4-86) 
a 



y/(l ~ Xl ) (1 - X 2 ) (1 - X 3 ) (1 - x 4 ) = (C 1)( " 2 , 

v /(C-x 1 )(C-x 2 )(C-x,)(C-x i ) = - (C ~ 1KC " 3 . 



(4.87) 
(4 



which immediately simplifies the two-cut spectral curve as 

v(z) = - z(z-m*-a (4 ' 89) 

Further notice, as we have discussed before, that above we have 3 equations for 4 
unknowns, and there is still one further moduli to consider; either t\ or t 2 (they are 
not independent as h +t 2 — t), 



1 1 3 

a+ -Jjoi, (4.90) 



3 

2 " ' 2 

1 1 
- a H — 

2 2 

8=0 



i=Q 

3 



J2 «i- ( 4 - 91 ) 



4.5 Some perturbative calculations 

We shall now discuss some perturbative calculations in the matrix models considered 
above. By perturbative, we mean a 't Hooft expansion in g s . As a warm up we focus 
on the matrix model (I4.44p . One can compute the partition function exactly (see 
appendix |A.2[) : the result can be written 

= G 2 (N + 1) G 2 (-N - 2a /g 8 + 1) G 2 (N + 2 Ql /g 8 + 1) G 2 (N + 2a 2 /g s + l) 
(2tt) n G 2 (-2a /g s + l) G 2 (2a 1 /g s + l) G 2 (2a 2 /g s + l) 

(4.92) 
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where G 2 (z) is the Barnes function. The Barnes function has the asymptotic ex- 
pansion 

logG 2 (iV + l) = ^N 2 log N+^N log 2vr-^iV 2 -^ log N + ('(-l) + 

+0 ° 1 B 
+ ^N 2 ^ 2 2 g{2g -2)' (4 ' 93) 

with 5 2g the Bernoulli numbers. However, (14.931) only deals with the N — > +oo, 
or g s + , asymptotic region. Depending on the sign of the finite parameters a\ 
and a 2 , as g s —> + one will have — — either going to +oo or to — oo, and one thus 
needs to also understand the asymptotics of the logarithm of the Barnes function 
in the region N — > — oo. However, it turns out that, from the point of view of the 
perturbative expansion, this sign difference is not very relevant. To clarify this issue, 
first notice the relation 

logG 2 (l-iV) = logG 2 (l + iV)-iVlog27r+ / da; ixx cot ttx. (4.94) 

Jo 

Explicitly evaluating the integral we find 

17T 

logG 2 (l-iV) = logG 2 (l + iV)-iVlog27r + — (1-6A^ 2 ) + 

+iVlog (1 - e 2 ^) - ^- Li 2 (e 2 ^) , (4.95) 

where Li 2 (^) is the dilogarithm. The logarithmic and dilogarithmic contributions 
can be expanded as 



N log (1 - e 2 * iN ) - — Li 2 (e 2 * iN ) = - V (- - 

&v ; 2vr v ; ^ \m 2mm 2 ) 



e 27 ™. (4.96) 



As explained in [32], this is actually the instanton contribution to the Barnes func- 
tion, also describable in terms of Stokes phenomena (across the +| Stokes line). In 
other words, this contribution is purely nonperturbative and we shall neglect it at 
this stage, i.e., from a purely perturbative point of view we may use following result 

17T 

log G 2 (1 - N) ~ log G 2 (1 + N) - N log 2tt + — (l - 6N 2 ) . (4.97) 

From the above discussion it follows that the logarithm of Z, in the 't Hooft 
limit, has the expansion^ 2 ! 



F = \ogZ = Y J 9l 9 ~ 2 F 9 , (4.98) 



12 As usual, we are only considering the real part of the free energy. 
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where, if we define T = F — ^N 2 \ogg s — |iVlog27r + logG 2 (iV + 1) (essentially 
amounting to the Gaussian normalization of the free energy) and use the relation 
(jOgD with k = 2, we find 

Jo 3pf = ^ (-oo + «i + a 2 ) 2 (\og[-a + «i + a 2 ] - - 2a 2 , ^log[2a ] - ^ 
+ i (a - an + a 2 ) 2 Uog[a ~ a x + a 2 ] - ^ J - 2a\ Hog[2ai] - ^ 



1 n2 / 1 „ r „. , „. 1 3 ^ o„2 A r_ , 3 



+ - (a Q + ai - a 2 ) ^log[a + a x - a 2 ] - - J - 2a 2 ^log[2a 2 ] - - J (4.99) 
for t7 = 0; 

3 f 1 — a + a x + a 2 1 . a - ai + a 2 1 a + a x - a 2 

>i = log log log 

1 12 s 2a 12 s 2ai 12 s 2a 2 

(4.100) 

for g = 1; and 



x-3 P f 

9 " 2^ (2^7 - 2 



-Qf + ai + a 2 ) 2 " 29 - (2a ) 2-29 (4.101) 
+ (a - ai + a 2 ) - (2aiJ + (a + a x - a 2 ) y - (2a 2 ) 



for g > 2. As alluded to above, it is also rather straightforward to compute the 
full nonperturbative contribution to this result. Because this is far from our present 
discussion we refer the reader to (12] for details, but the main idea essentially follows 
from the application of 

disc log G 2 (N + 1) = i g (M + _M e -*r|*|m (4.102) 

to the expression for the free energy (where the discontinuity of the free energy will 
yield the full tower of multi-instanton corrections). In this case one simply obtains 

^. j. i / 27T (— «0 + Oil + a 2 ) i^g \ 2ir(-a + a 1 +a 2 )m 

3pf 27r ^£^v m ™ 2 / 



m=l 

+00 



1 / 2lT (a — ai + a 2 ) # s \ Map-gi+aajm 

2 ^^V m m V 



+00 



27r («q + on ~ a 2 ) g s \ 2n( ao + ai - a2 ) 



- t^-XJ — : — ; — " • -1 



m=l x 



2 

+00 / , _ \ . +00 



47ra g s \ ^apm i ^-^ / Arcai g s 



i / 47ra 2 | g s \ e _*sgas (4 103) 
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Notice that this is an exact result, to all loops and including all instanton numbers. 

Even though we have (for this case) an exact perturbative expression it is useful 
to also compute Fi using a method that generalises to more complicated situations 
where exact results are unavailable. There is a universal formula for F% which takes 
the form [4Tpl 

F x = -llog (M(a)M(b)(a - b) 4 ) (4.104) 
where (this follows from ( I4.30P and (14.54)) ) 



a 



M(a)M(b) = 4^^ • ( 4 - 105 ) 



It immediately follows, using (14.54)) . that 

Fi = -— log(2a + 2ai + 2a 2 ) - — log 

"T^ log ~ n T^ log 7, ' ( 4 - 106) 

12 (X\ IZ a 2 

reproducing the result we have previously obtained, up to some irrelevant numerical 
terms (this expression explicitly includes the Gaussian contribution — y^logt). 

Next we turn to the case corresponding to the chiral four-point function in the 
Liouville theory, namely the matrix model with potential (14. 70 p . In this case we do 
not have an exact solution, but as discussed in the section 14.41 we can get tractable 
expressions if we work order by order in (. 

As an example we consider the special case discussed at the end of section 14.41 
and focus on F±, which may be computed in a straightforward fashion from the 
universal result (I4.104p . We have 

, w \ 2cti 2a 2 2a 3 

M ( z ) = ~^= + l x 77 w + 



(z - 1) y/(l-a)(l-b) (z - C) y/(C-a)(C-b) 
2«i + 2a 3 _^ 2a2 ^ 2a 3 ( + 2a 3 (a + b) ( + ^ ^ 



zy/ab {z - 1) y/(l - a)(l - b) z 2 y/ab 2z (ab) 
2a ° f l+ / 3C , +0{( 2 )), (4.107) 



z(l - z) \ z (ai + 0:3) 



leading to 



M M M M = m 4 1i m f 1 + ^ + (« 2 )) • < 4108 > 

ab(l — a)(l — b) \ ab a\ + a 3 v / 



and further computing 



a 



fl-^^ + 0(C 2 )V (4.109) 



13 The universal expression was derived for polynomial potentials, but our results indicate that 
it also holds for the multi-Penner-type potentials considered in [T]. 
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where Q was defined in (I4.83p . and 

b + a = al + {ai + a3)2 ~ al + al ~ ( 2 a ; + a3)2 ± « 3 c + o (c 2 ) (4.1 io) 

it finally follows, after putting all the above expressions together, 

M(a)M(b) = Aa [ ^ U + ^^ + 0(e)), (4.111) 

(f"«) 4 = T fl " J^T- + O (C 2 )) , (4.112) 
al \ 2ai + 2a 3 'J 

and (here we also included the ( 2 terms) 

1 , 2Q , . 

= log (4.113) 
12 a ("l + "3) "2 

aia 3 (— a +«i— a2+«3)(«o+«i— tt2+«3)(— ao+«i+«2+«3)(«o+«i+«2+a3)^2 

32( ai + a 3 ) 6 

Notice that, up to numerical constants that we drop, and to first non-trivial order 
in £, this result is exactly the same as the result in (I4.106p . except for the shift 
aci — > aci + «3. Note also that the order 0(C) contribution to F\ vanishes. 

The above result can be compared to the corresponding result for the Nekrasov 
instanton partition function of the SU(2) theory with four fundamental hypermul- 
tiplets. In this case one easily obtains 

F inst = -aVi(m 2 ) + 2a 4 a 2 (m 2 ) - 3aV 3 K) + 4tx 4 (m 2 ) ^ 2 + 3 

128a^ 

where a^{m 2 ) = Y^h<-<i m ?i ' ' ' m 'i k - Note that the 0(y) term vanishes in agree- 
ment with (14.1131) identifying y with (. Furthermore, after implementing the rela- 
tions, 

rrii = di + «3 , ni2 = oli — «3 , m 3 = «o + «2 , fn^ = «o — ol% (4.115) 

together with a = rrii we see that also the second order terms in (14.1131) and (14. 114[) 
agree perfectly. (The definition of a in the Nekrasov expression differs from the one 
in (14.731) .) This result is consistent with the analysis in section I4T21 and supports the 
approach in pp. 

Computing JF g for g > 2 would require heavier machinery, see e.g. [15] , and will 
not be attempted here. Similarly, an explicit example involving the full-fledged two- 
cut geometry would take us too far afield. The integrals of the periods associated to 
the 't Hooft moduli are generically hard to evaluate exactly (although it is possible 
to do so in the present situation) and even harder to invert in order to find explicit 
solutions for the endpoints of the two cuts. A possible way out is to resort to 
perturbation theory, along the lines of [371 EE] , but we shall leave this question for 
future work. 
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Let us close with a comment about the relation to the chiral four-point function 
in the Liouville theory to leading order in (. In order to reconstruct the full Liouville 
three-point function at this order pertubative calculations are not enough, one would 
also need to add the full set of nonperturbative corrections and in general also let 
1, in order to obtain the desired result. It seems plausible that a matrix model 
perturbation theory in ( exists (for general (3) which is exact in g s order-by-order, 
but we leave this question for future work. 



5 The A r matrix models 

In this section we perform several calculations in the A r quiver matrix models. The 
resulting expressions are compared to the corresponding expressions in the A r Toda 
theories and the A r quiver gauge theories. 



5.1 The three— point function 

The above analysis of the A\ matrix model three-point function (see section 14. ip 
can be extended to the A r theory for any r. The relevant integral is 

H (i,I)<(j,J) iJ h 1 

Here ot\ = (ax, e») with ax = a\Ai and similarly for a 2 (see appendix lA.ll for our Lie 
algebra conventions). The integral (15 .ip can be explicitly evaluated provided one 
imposes the restriction a\ = x5 ir (i.e. ax = xA r ); the result is 



JT^JJ " r([ai + ---ai]/ex+j-i + l + (I-l+i-m 



l<i<j<r 



r([c4 + • • • a> 2 + a[}/£x + j - % + 2 + (/ - 2 + i - j + N, - N j+l )f3) 

_ jj -q emm + i + Ni+i - mnm (52) 
»=i i=i 

where A^ = A^ r+ i = 0. The result (15.21) depends on a very particular choice of 
integration contour, which is quite subtle and will not be discussed here; see jH] for 
further details. 

Using various Lie algebra results (see appendix IA.1I for further details) one can 
show that the relation 

a\ + a\ - al - £2A ij N j = , (5.3) 
where ax = xh r (i.e. a\ = >c5i r ) implies that 

- e 2 (N i - N l+1 ) = -Ar X K + A^ l r x + (a 2 , h i+l ) - (a Q , h i+1 ) . (5.4) 

(Note the special case, e%N\ = + (a 2 , hi) — (a , hx)-) Furthermore, 

a\ - A~ l x + A-^h = -Jy . (5.5) 
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Using these results together with 

a\ H h ol\ = -(a 2 ,Uj+i ~ Ui) , % - j - 1 = ~(p,Ui - u j+ i) , (5.6) 

and ( 14.41) . it follows that the expression ( 15.2ft can be written (we suppress an unim- 
portant prefactor and T 2 (x) is short-hand for T 2 (x| — e±, —£2)) 

-j-r T 2 ({a 2 +pe, hj}-{a +pe, h i )+^ l )r 2 {-e-(a2+pe, hi)+(a +pe, hj}-^) 
i<i<j< r +i r 2 (-(a 2 + pe, e ii ))r 2 (-e - («o + pe, ey)) 

r 2 (-e-(a 2 +pe,fe 1 ) + (a +pe,fe 1 )-^ I ) r ^ x 

x =n; \ T 2 ((a 2 +pe,h i )-(ao+pe,h i }^ — — ) 

r 2 (-£-x) r+1 

n r 2 (- £ - iVj£ 2 ) • (5 ' 7) 

The first thing to note is that the factors on the last line is a phase. As in the A\ case 
discussed in section 14.11 the above expression can be compared to the three-point 
function in the A r Toda theory. After suitably rescaling the vertex operators with 
multiplicative factors depending on their momenta, the A r Toda theory three-point 
function (12.9ft can be written 

r+l 



X 



X 



X 



II ^ + {ai- Qp, h) + {a 2 -Qp,h j })\ 

l<i<j'<r+l 

II \ T b{Q - — y - («i - Qp, hj) - (« 2 - Qp, hi}) \ 2 

l<i<j<r+l 

r 6 (Q — - (ai - Qp, /ii) - (« 2 - Qp, M) (5i 

'r + l 

■r+l 



where r&(a;) is short-hand for T 2 (x\b,l/b). By similarly rescaling of the matrix model 
vertex operators we are left with the factors in the numerators on the first two lines 
of (15. 7p . After using the identification ( I4.8p together with Qp — a — > —(Qp — «i) 
we find complete agreement with the "square-root" of (15. 8ft . The expression (15. 7ft 
should probably also be related to (the perturbative piece in) the T r+ i = 73,0 (A-) 
theory where one of the masses satisfies a restriction inherited from a% = xA r via 
the AGT relation. 

We note that in [251 ES] a complex version of the above integral was used to 
derive the three-point function in the A r Toda theory when one of the momenta 
takes the special value xA r . 

In the above evaluation of the integral (15.11) the condition a\ = xA r was imposed 
with a 2 left arbitrary. In (35] the above integral was evaluated for the rank 2 case 
with the alternative restriction: a 2 = xA 1 — (x + e)A 2 , with a.\ left arbitrary. 
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To understand why there is more than one possible choice which allows for an 
explicit evaluation of the above integral, we recall that in the A 2 Toda theory the 
condition a = xA r translates into the fact that the corresponding W primary state 
satisfies (see e.g. [26j [5] and references therein) 



3w(a) 
2A(a) 



L. 



\a) = 0. 



which implies 



A(a) 
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22 + 5c 



A(a) + - 



-w{a) 



where 



and 



A(a) 



(2Qp — a, a) 



w\a) 



18 



(a - Qp, h x ) (a - Qp, h 2 ) {a - Qp, h 3 



(5.9) 



(5.10) 



(5.11) 



(5.12) 



22 + 5c 

In (I5.12p the hi are the weights of the fundamental representation of the A 2 Lie 
algebra, cf. ( 1A.2I) . If we write a = a 1 Ai + a 2 A 2 , then it turns out that there are six 
one parameter solutions to (15.101) : 



a = xAi , a = xAi + 2QA 2 , a 
a = xAi — (x — Q)A 2 , a = xAi — (x 



kK 2 , a 
3Q)A 2 



2QAi + xA 2 



(5.13) 



Using that in our conventions e = —Q, we see that both the conditions (12. 8ft as well 
as the condition used in [15] belong to the set C I5 . 131) . In addition to these three 
solutions there are three more which differ from the other ones only when Q / 0. 
It is known that in the A 2 Toda theory these six possibilities do not correspond to 
distinct states, rather they are related via the so called shifted Weyl group acting 
on the momenta which changes the corresponding vertex operators by the so called 
reflection amplitudes^). Therefore there should also be a simple relation between 
the corresponding matrix integrals; in particular, if one can be explicitly evaluated 
then that should also be the case for the others. 



5.2 Higher— point functions 

As in the A\ case we can analyse a certain class of four-point correlation functions 
exactly. We need the following result 



' N A TT FT [a\ + --- + a\]/e l + i + (N x - I - i + 1)0 



fN \ 



i /=i t a 2+ •••+«! + a\]/£ + i + l+(N 1 +N i -N i+1 -I-i)f3'- 



14 We thank Yuji Tachikawa for clarifying this point. 
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where the A r Selberg integral S r (ai,a 2 , (3) was defined in f l5.ll) . a\ = >c5 ir , and 
e^(Ai) is the £th elementary symmetric polynomial defined as 

e/(Ai)= Af---A?. (5.15) 

h<-<le 

From the result (15.141) it follows that 

S P (a 1 ,a 2) /3;z)= /jl^ J] 1^ ~ A/| M « JJ(Af)^(l - Af)^T JJ(^ - A{) 

il (i,I)<(j,J) i,I I 



r Ni-i 

X 



-Q -Q [a 1 2 + --- + a i 2 \/e 1 + i + (N 1 -I-i + l)f3 



1 [«2 + • • • + «2 + + * + 1 + (^1 + N i ~ N i+1 -I-i)P 

Nl 7 l N I 

a 2 , /?; 0) £ ^(-if^Tyy (5-16) 



£=0 

JL £ll -£(a£ H h aj + «1 + iei + ei) + 1 + iV< - N i+1 - i 

x H!! -H a i + • • • + «2 + ^i) + j - * + 1 ' 

As above one can show using various Lie algebra results (see appendix I A . 1 1 for further 
details) that 

-£ 2 (Ni-N i+1 ) = -A- 1 s 1 +A^ hl e 1 -A; r 1 x+A^ hr x+(a 2 ,h l+1 }-(a ,h l+1 } . (5.17) 
where we have used a± = xA r and «3 = SiAi. Furthermore, 

- A 7 ! 1 + A 7+i,i = , <*[ - K 1 * + A 7+i,r* = ■ (5-18) 

Using these results together with 

oc\ H hag = -(a 2 ,/ii+i - /ii) , i = -(p,h i+ i - fit) , (5.19) 

it follows that 

6 r (ai,a 2 ,p;2) = S r (ai, a 2 , p; 0) 2^ ^[ (B 1 ) e ---(B) e 

= S r (a 1 ,a 2 ,0;O)r+xFM 1 ,...,A r+1 ;B 1 ,...,B r ;z) (5.20) 
where (X) n = A(X + 1) • • ■ (X + n — 1) is the Pochhammer symbol, 

B % = — -(a\ H h a\ + iei) - i + 1 = — (a 2 + ep, h l+l -hx) + l (5.21) 

£2 £2 
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and 



A = ( — + £ i— 7T + ( a 2 + ep, h x ) - (a + ep, hi) ) . (5.22) 

Note that (-I)^t^W = (Ai) e and that the sum over £ in (I5.20p can be extended 
to oo since (Ax)i — for I > Ni. 

To compare the result ( 15.20P with the corresponding result in the A r Toda theory 
we recall that in [251 126] it was shown that the correlator 

(V-bAi (z)V ai (0)V a2 (oo)T4 Ar (l)) (5.23) 

satisfies a differential equation of hypergeometric type whose solutions involves 
r+1 F r (Ai, A r+ i; Bi, . . . , B r ; z) where 

Ai = b( — ^— -b—^— + (ai-Qp,hi) + (a 2 -Qp,hi)] 
\r + 1 r + 1 / 

B t = l + b(ai-Qp,hi-h i+ i) (5.24) 

Replacing ati — > a 2 , a 2 — Qp —> — (« — Qp) and using the rule fl4.8[) we see that 
fl5T2"4"j) agrees perfectly with f l5T2T|) and (ET22)) . 

One can also show that the Nekrasov partition function leads to the same result, 
see [7| for a discussion. It should also be possible to analyse the correlation functions 
involving the insertion of V_Ai/6 as we U as the case with several insertions of VL&Ar 

To analyse general correlation functions without imposing the above restrictions 
on the momenta is much more difficult. However, we stress that matrix model per- 
turbation theory can handle any correlation function, although this method appears 
to be somewhat cumbersome and it is not clear what the meaning of the resulting 
expansion is in the 2d CFT. 

5.3 The curve 

We now turn to the discussion of the loop equations and the large N matrix model 
curve. In pQ and in sections 14.31 and H~4l some examples of curves in the A\ case were 
presented. Here we mainly focus on the A 2 case. We call the two matrices $ and $ 
and the associated potentials W and W. The (non-hyperelliptic) curve is known to 
be of the form [TBI EI] 

x 3 = r(z)x + s(z) , (5.25) 



where 



r(z) = -\W\zf + W\zf + W\z)W\z)} (5.26) 
3 

_ s , / tr (W(z)-W'(*)\ \ _ gt L ( *"« - «"(*)' 



2-$ / / \\ Z - $ 
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and 

s(z) 



[W'(z) + 2W'(z)}[2W'(z) + W'(z)][W'(z) - W'(z)\ 




W'(z) - W'($) 



z — $ 
W'(z) - W'(§) 



z 



+9s < tr 



where 



UJ r (Z) 



: (2W\z) +W\z)) 



uJz) 



: {W\z) + 2W\z)). 



(5.27) 



(5.28) 



(The expressions on the last two lines in (15.27ft can be simplified in the eigenvalue 
basis by using the saddle-point equations, but we will not need the resulting expres- 
sion here.) 

Inserting the explicit expressions for W and W as sums of logarithms^ 
v v 
W($) = ^m i log(z i -$), W{$) = ^m i log(z i -<£), (5.29) 



i=i t=i 
and using the above expressions for r(z) and s(z) we find 



x 



3(2) 
2p-2 



X + 



p£Uz) 



nu(z-zi 



(5.30) 



where P^p^iz) and P 3p i 3 (^) are polynomials of degree 2p— 2 and 3p— 3, respectively. 
This curve is of precisely the right form to agree with the expression in jl] (which 
was obtained by starting from the earlier result [47J). However, there is a further 
property that should to be checked. 

Recall that, in the A 2 case, there are two types of punctures, one full and one 
basic. For the basic (or special) puncture there is a relation between m, and fhi. In 
our conventions that relation is fhi = 0. Now for a special puncture at z = Z{ there 
is a relation between P^(z) and Pw(z) that has to be satisfied at that location, 
see (3.25) in [4J. To check that this condition holds for the matrix model curve it is 
sufficient to focus on a single special puncture which we can take to be located at 
Zi = 0, i.e., W($) = mlog($) and W($) = 0. We need to check that 4r (z) 3 -27 's(z) 2 
scales like -4 (naively it would scale like -§■). In other words we need 



A{-[W\zf + W'{z) 2 + W\z)W'(z)]f 
3 

W(-tL\W'(z) + 2W\z)}[2W'(z) + W\z)][W\z) - W\z)]f 



15 Since the curve is derived in the limit of vanishing e the AGT relation is ai 
suitable definition of the masses. 



(5.31) 



wii using a 
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and 



12(-[W'(zf + W'{zf + W'(z)W'(z)}) 2 (f(z) + /») (5.32) 
+ 2[M/X^)+2^ , (^)][2M/'(^)+^ / (^)][M/X^)-^ , (^)]K(^)/>) -u> r (z)f(z)) = 0, 
where 

(5.33) 

Both the above equations are easily shown to hold for the special puncture with 
rhi = thereby establishing the equivalence with the results in [4"ffi 

For higher rank curves the general structure of the matrix model curve is also 
known; see [48J for the state-of-the-art knowledge. The matrix model curves can 
be analyzed and compared to the gauge theory curve as above. 



6 5d gauge theories and q— deformed matrix models 

Nekrasov partition functions can also be defined for supersymmetric gauge theories 
in five dimensions formulated on R 4 x S 1 [3J. As an example, in the SU(2) theory 
with four matter hypermultiplets in the fundamental representation the instanton 
partition function is 

7 _ \ ^ \Y\ TT TT 'Pjo-mi Ym, s)) / R \ 

Anst- 2^V 11 11 g( * _ a y Y s)(£(a -a Y Y s) — e) 1 J 

where (0,1,0,2) = {a, —a) and 

£(x,Y m ,Y n ,s) = 8mh(R(x-e 1 L Yn (s) + e 2 (A Ym (s) + l))), 
4 

V{a m ,Yi,s) = Y[smh(R(a m - (j - l)e 1 - (i - l)e 2 - m f )) (6.2) 
/=i 

with R the radius of the S l . See section [2T21 for more details about the notation. 
The expression ( 16. II) can also be obtained from topological string considerations see 
e.g. [32] [50]). The partition function also has a perturbative piece whose explicit 
expression can be found e.g. in [3l I5T]. 

The question we would like to address in this section is: Are there matrix model 
and CFT descriptions of partition functions of the type (16. ip ? 

In the recent paper [TTJ a proposal was made for the CFT description of the pure 
SU(2) theory in five dimensions. This proposal involved the so called g-deformed 
Virasoro algebra [52] ■ It is very natural to expect that there is a CFT description of 



16 There is also another solution to (|5.31[) and (|5.32l) . viz. — rhi. This solution is precisely the 
alternative solution discussed at the end of section IQl (note that the curve is derived for e = 0). 
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5d (conformal) quiver gauge theories which involves g-deformed Virasoro [52] arid q- 
deformed W algebras [53] . Unfortunately the representation theory of these algebras 
is not very well developed. Analogues of the primary fields and their quantum 
numbers are known, but the analogue of e.g. the relation 

[L m ,V a ] = z m [(m + l)A(a)V a + z(L^V a )} 

= z m (mA(a)V a + [L ,V a }), (6.3) 

is not known (as far as we know). This fact complicates the analysis and makes 
direct calculations of chiral blocks difficult. 

Instead, we try to obtain a matrix model description. A natural starting point 
is to look for a generalisation of the Ai three-point function (14. ip . 

There is a known g-deformation of (14. ip in the literature, which can be written 

where < q < 1 and J dx q is the so called g-integral (or Jackson integral) defined 
via 

/ d q xf(x) = (l-q)J2f(q k h k - (6-5) 

In the limit q — > 1 _ this expression converges to the Riemann integral f(x)dx. 
Furthermore, (a; g)oo = YYkLoi^- ~ a 1 k )y an d we a l so use the notation: 

(a; q) e = (1 - a)(l - aq) ■ - ■ (1 - ag^ 1 ) = 7 ^f- . (6.6) 

Based on the above expression we tentatively propose the rule 



In the limit q — >• 1 , V£(z) — > V a (z). In the special case z = we assume that 



V*(0) = V a (Q) = n / (A / ) 2a/£l - Note that when a = ej2, V«(z) = Ui z TTp^ 3h 



reduces to Yli( z ~ g _1//2 A / ). 

The above three-point function can be evaluated (at least when /3 is an integer) 
and leads to a product of g-gamma functions, but since we are not aware of a 
g-analogue of ( I4.4|) we will not discuss the result here. 

Instead we turn to the four-point function with 0:3 = £i/2. Using the result in 
the resulting expression can be explicitly evaluated 

S"(a l ,a 2 ,(3;z)= / [[d q X ( z -—)(\ I )n — ||(A ) ^ — ^- 

J 1=1 q (g £1 A/;g)oo kj KTxj)™ 

= S"(ai, a 2 , /3; 0) 2 1 (?- /3iV , ^Ai+^Ai+^W-^. ^ ~ z) (68) 
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_ 2a 2 /e 1 -l 

where z = zq 2 and 



a ! a u n \ (A; q)k{B\ q)k k n \ 

MA,B-,C-, q ,z) = ^ {c ^ q)k Z . (6.9) 

(Actually, since iV is an integer 2 0i in the above expression reduces to a so-called 
little g-Jacobi polynomial.) 

The integral (16.81) was also discussed in [55], albeit in a somewhat different guise. 
In that paper it was shown that the integral satisfies a certain difference equation. 
This is a g-analogue of the result in [36] (c/. the discussion in section H~2l) . However, 
the analysis in [M] is more transparent, although the case corresponding to multiple 



insertions of V£, 2 is not discussed in [5 
The above result (16. 8p can be rewritten 

2 MQ- P ", g P«MW^-ffl ;g C*«Mi) ; gP } q2 J £i _ 1/2 ) (6.10) 

_ " tl sinh(R[-N + k]) sinh^- 1 ^ + jf + 2) + iV - 1 + k}) 
~ ^ sinh( J R[/3- 1 (2ai/ei + 1) + fc]) sinh( J R[l + jfe]) {ZQ ' ' 

where we have used g = e~ 2R ^. 

To compare this expression with the one arising from the Nekrasov partition 
function, we go through the same steps as in section 14.21 and make the choices 
mi = cl and vri2 = —a + e\ which implies that (16. ip reduces to 

^ t A. s inh(i?[(a 1 +a 2 +a 4 -^-])/e2 + k) sinh(i?[(ai+o;2+a4-|ei-e2)/e2 + k\) 
U V \} sinh( J R[( ei -2a 1 )/e 2 + fc])sinh( J R[l + A;]) ' 

(6.11) 

where we have also used the AGT relation (14. 16ft . This expression is readily seen to 
agree with (I6.10p using y = zq 1-13 and the same arguments as in section (14. 2p . cf. 
the discussion after eq. (I4.18p . 

In [55] there is also an extension of the above result to the case with multiple 
insertions of Vff, T In this case the result involves the function 2 < £ ) i 9 ' t ' ) (A B; C; z) 
which is a q analogue of ( I4.2ip and involves Macdonald polynomials rather than 
Jack polynomials, see [55] for further details. So far we have only considered the A\ 
case; it should also be possible to consider the A r case using the results in |44j. 

We close this section with a few words of caution. There are in general several 
possible g-deformations and the one above may not be the right one. Also, we 
should mention that in [51] another deformation of the matrix model was shown to 
be related to Nekrasov partition functions for five-dimensional gauge theories. The 
deformation in [56J replaces the Vandermonde determinant with ri/<j sinh(A / — A J ). 
This possibility was also mentioned in [1]. 

7 Discussion and outlook 

In this paper we have studied the A r quiver matrix models which were introduced 
in [Tj and argued to capture correlation functions (chiral blocks) in the 2d A r Toda 
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field theories and Nekrasov partition functions (instanton partition functions) in the 
4d A r quiver gauge theories. 

From the point of view of the matrix model the expansion in Zi (the locations 
of the vertex operators) is somewhat awkward, but we have shown that several 
known results can be rederived from the matrix model; some of our checks are 
quite non-trivial. It would be interesting to develop the matrix model technology 
further, and, for instance, to clarify the choice of integration contour and to develop 
a perturbation theory in z%. 

We also made a proposal for an extension of the matrix model to capture the 
Nekrasov partition function of 5d quiver gauge theories. This speculative proposal 
passed a non-trivial check, but deserves further study. 

One open problem is to extend the analysis to the other ADE Lie algebras. 
Let us make a comment about the D r case. The matrix model curve for the D r 
model can be extracted from |38], eq. (3.51) and can be seen (after some changes of 
notation) to be of the same general form as the curves in [57] . It would be interesting 
to study this in more detail. 
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A Appendix 

A.l A r roots and weights 

Here we collect some standard results for the A r Lie algebras. The root /weight space 
of the A r Lie algebra can viewed as a r-dimensional subspace of W +l . The unit 
vectors of IR' r+1 will be denoted «j (i = 1, . . . , r + 1) and satisfy (ttj, Uj) = 5ij. The 
simple roots are e, = Uj — u i+1 (i — 1, . . . , r) and the positive roots are = «j — Uj 
(with 1 < i < j < r + 1). The Cartan matrix is A^ = (ei,ej), and its inverse is 
A# = 7+1 min(z, j) [r + 1 — max(i, j)]. The Weyl vector, p, is half the sum of the 
positive roots; hence p = ~ Y7i=i( r ~ 2z + 2)u{. The fundamental weights, Aj, are 
defined as 

r+1 

A 4 = « 1 H hM *~I7T rX^'' ( i = 1 T--' r ) (A- 1 ) 

r 3=1 

and satisfy (A^ej) = dij. Note that Y^i=\ ^-t = P- Finally the weights of the 
fundamental representation can be chosen as 

1 ^ l 
hi=Ui- -—^^Uj = A 1 -^2e j , (i = 1, . . . ,r + 1) (A.2) 

j i =1 
Note that hi = Ai and ^ . hj = 0. 
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A. 2 Orthogonal polynomials and the A\ three— point function 

Here we present an alternative evaluation method for the integral (14. ip in the case 
(3 = 1 using orthogonal polynomials. Consider the one-matrix model partition 
function 

„ N 

z = Tspjff / n n>' - ^ •* E£i ww • m 

v 7 ' J 1=1 KJ 

and introduce orthogonal polynomials, {p n (z)}, with respect to the measure 

dfi(z) = e^ w{z) dz (A.4) 

as 

dfx(z) p n {z)p m (z) = h n S nm , n > , (A. 5) 



where one further normalizes p n (z) such that p n (z) = z n + ■ • • . Noticing that 
the Vandermonde determinant H/<j(^ / ~~ ^ J ) 2 equals detpj_i(A 7 ), the one-matrix 
model partition function may be computed as 

V ; n=0 

In the case of interest to us, the potential is 

k 

W{z) = tr 2a * - z ) > ( A - 7 ) 

a=l 

which, in principle, forbids the use of standard orthogonal polynomial techniques. 
However, the fact that the non-polynomial structure is logarithmic actually allows 
us to get around this issue when k = 2, as we shall see now. Indeed, in this case 
(setting zi = and z 2 = 1) the measure associated with (1A.7|) becomes 

d(i(z) = (1 - z) 2a2 z 2a ' dz } (A.8) 

and is immediately related to the orthogonal polynomial family of Jacobi polyno- 
mials. The combination 



j^) {z) s n!r(n + q + 7 + l) ^ _ 
n K J r(2n + a + 7 + l) n V ; 



(A.9) 



where P„ (#) is a Jacobi polynomial, is normalized such that Jn (z) = z n + • • • 
and satisfies ^ 

[ dz(l-z) a z\J^\z)J^\z) = h n S nm , (A.10) 
Jo 

with 

_ T (w + a + 1) T (n + 7 + 1) T (w + a + 7 + 1) 
n n ' r(2n + q + 7 + 2)T(2n + q + 7 + l) ' 1 J 
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Using (1A.6|) this immediately leads to the exact result 



j-j^T (n + 2a 2 /g s + 1) T (n + ai/g s + 1) T (n + a 2 + ai + 1) 



(2^ ^ = A ' r (2n + a 2 + «i + 2) T (2n + a 2 + «i + 1) 



1 TT»' r(n + a2 + 1)r(n + ai + 1) (A19) 
(2vr)^ 11 ' r(iy + n + l + « 2 + ai ) ' 1 ' 1 

v ' n=0 

which agrees with (14.21) (when (3 = 1). 
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